A freely relcixing polymer remembers how it was straightened 
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The relaxation of initially straight semiflexible polymers has been discussed mainly with respect 
to the longest relaxation time. The biologically relevant non-equilibrium dynamics on shorter times 
is comparatively poorly understood, partly because "initially straight" can be realized in manifold 
ways. Combining Brownian dynamics simulations and systematic theory, we demonstrate how 
different experimental preparations give rise to specific short-time and universal long-time dynamics. 
We also discuss boundary effects and the onset of the stretch-coil transition. 
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The mechanical properties of semiflexible polymers, 
which form an integral part of the cell structure, arc 
of high relevance to the understanding of cell elasticity 
and motility [l], 0| ■ External stress not only remarkably 
changes static and dynamic features [1, 0, 0], but has 
also important biological implications, e.g., for enzyme 
activity on DNA 0, [l]. Particularly intruiging as- 
pects of stress-controlled behavior can be observed for the 
relaxation of semiflexible filaments from initially nearly 
straight (i.e., highly stressed) conformations. In recent 
years, many experimental and theoretical studies have 
addressed this paradigmat ic p roblem of polymer rheol- 

M,M, M, H, m, iMlllllllalf)^ often primarily fo- 
cused on the influence of hydrodynamic interactions on 
the longest relaxation time tR, which is a key identifler 
of the stretch-coil transition. However, on much shorter 
times the polymer dynamics is predominantly controlled 
by the hi ghly nontrivial internal conformational relax- 
ation [2^. l30l|. which plays a relevant role in many bio- 
logical situations ranging from the viscoelastic response 
of polymer networks 31 1 to molecular motor kinetics (3^ 
and DNA supercoiling dynamics [13, HI] . This aspect of 
the relaxation is still poorly understood, the more so as 
standard analytical techniques based on linearized equa- 
tions of motion fail due to inherent nonlinearities initi- 
ated by strong perturbations [H, [sl]. Further, because 
a completely straightened polymer conformation can in 
practice not be realized in the presence of thermal noise 
from the environment, the short-time dynamics of an ini- 
tially "nearly" straight filament will reflect the way it was 
straightened: filaments can be stretched by optical tweez- 
ers [2^ [2^, by electric fields [isl . [TtI . [l^. .23.. ,36il. o r by 
flows of different geometry d |l3 [ll [la [!£ [Illli, [sj, 
but a straightened contour can also result from low ini- 
tial temperatures. In any case the relaxation dynam- 
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ics is driven exclusively by stochastic forces. This raises 
the question how results obtained with different setups 
should be compared and when the dependence on initial 
conditions fades out. 

In the following, we present results from computer sim- 
ulations combined with a thorough and exhaustive the- 
oretical analysis to explain how fundamental differences 
in the short-time relaxation emerge from different ex- 
perimental preparation methods but give way to univer- 
sal long-time relaxation. Four idealized initial conditions 
(see the cartoons in Fig. [T|) are shown to lead to quali- 
tatively distinct behavior despite superficial similarities. 
"Force" refers to mechanical stretching, i.e., a strong ex- 
ternal stretching force /pre that is suddenly removed on 
both ends, for instance in a setup using A-DNA, opti- 
cal tweezers, and restriction enzymes @- Secondly, the 
term "field" is used for experiments employing an elec- 
tric field [l3| of strength E for stretching, where one end 
is always kept fixed. Once switched off, such fields give 
rise to relaxation dynamics similar to the one in setups 
using homogeneous elongational flows Q of velocity v. 
Further, we denote by "shear" the stretching by planar 
extensional shear flows of shear rate 7 in a symmetric ge- 
ometry [HjH^I, see Fig. [TJc). Finally, "quench" refers to 
a scenario where the temperature is suddenly increased 
by a large factor 9 from a small value T/9 near zero 
to its flnal value T. This setup is more feasible for com- 
puter simulations, but the equivalent sudden drop in per- 
sistence length ip might be experimentally realizable by 
chemical reactions. 

The paper is organized as follows. In the flrst sec- 
tion, we present results from Brownian dynamics simu- 
lations for each of the four different setups. The next 
section presents a qualitative discussion of the underly- 
ing theoretical model resulting in scaling laws for perti- 
nent observables, which readily suggest intuitive explana- 
tions for the qualitative differences between the scenarios 
and their universal long-time asymptote. A detailed and 
somewhat technical derivation of these asymptotic scal- 
ing laws is contained in the third section, where we also 
analyze the effect of different boundary conditions. In 



FIG. 1: (color online) In the "force" scenario (a), a stretch- 
ing force /pre is suddenly removed at both ends. In the 
"field" setup (b), one end is held fixed and the electric field 
of strength E (or homogeneous elongational flow of velocity 
v) is switched off, similar to the "shear" case (c) with a sym- 
metric extensional shear flow of shear rate 7. In the "quench" 
experiment (d), the temperature is suddenly increased by a 
large factor 9 from an initially small value. 

the fourth section, we present a quantitative comparison 
between simulation results and theory. At the end of 
the paper, we discuss experimental implications includ- 
ing quantitative estimates of control parameters in typ- 
ical realizations, the onset of the stretch-coil transition 
and the influence of hydrodynamic interactions. 

I. SIMULATION RESULTS 

In the Brownian dynamics simulation, we employ the 
standard free-draining bead-spring algorithm for worm- 
like chains, were different environmental conditions dur- 
ing equilibration of the chains correspond to the four sce- 
narios introduced above. The equations of motion for a 
chain of total length L = Nb with + 1 beads of size b 
and mobility fi arc given by 

dtVi -Vi = -/iV,:[/ + r]i{t), (1) 

where the potential U = Ug+Uh+Uf contains a stretching 
part 

Ug = ^Y.i\^^^^-^^\-bf, (2) 

i 

a bending part 

i 

and an external potential U{. Here, ip is the persistence 
length, 7s is the stretching elastic constant, and U = 
is a normalized tangent vector. We use Gaus- 
sian noise with strength {r]i{t)r]j{t')) = Q^k^T 5i,j5{t—t') . 
The time step is 10~^to, where To = b^/{kBTfj,) is the 
self-diffusion time of the beads. The chains are equili- 
brated along the x-axis symmetrically to the origin un- 
der the respective stretching mechanism. In the "force" 
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FIG. 2: (color online) Quantitative results: change in pro- 
jected length (5i?||(t) and mean bulk tension fh{t) (inset) for 
computer simulations of a "force" (squares), "field" (trian- 
gles), "shear" (circles) and "quench" (diamonds) scenario, re- 
spectively. Simulation parameters were L = 2006, £p — 40b, 
7s = 6000/b, /pro = 10 keT/b, v = 0.18 b/ro, 7 = 0.0046/ro 
and = 35.7, such that -R||(0) ~ 0.971/ in all cases. 



case, Ui = — /pro(a^Ar — xq) and Vi = 0, while Uf = 
and Vi ~ (11,0,0)-'" for "field" setups. For "shear", we 
take Vi = j {xi, --yi, ^Zi)^ , and in order to prevent the 
polymer from diffusing out of the stagnation point, an 
additional harmonic potential Ui = \i iJ''^^ Xqqm drives 
the center-of-mass coordinate xcom back to the origin 
(cf. the feedback control system in Ref. [11]). In these 
scenarios, we equilibrate for lO'' time steps, while initial 
conformations are generated directly using the equilib- 
rium tangent correlations in the "quench" case. In all 
cases, J7f = and Vi = upon release. Ensemble aver- 
ages were taken over 150 realizations 

To characterize the relaxation dynamics, we concen- 
trate on two observables. One is the time dependent 
change 5R\\ (t) = R\\ (0) — R\\ (t) in the ensemble aver- 
age of the filament's end-to-end distance projected 
onto the initial longitudinal axis. Note that with this 
definition, 5R\\ is positive and increasing, while the ac- 
tual end-to-end distance shrinks during relaxation. The 
second observable is the mean tension ft,{t) in the fil- 
ament, proportional to the bulk stress a{t) in a poly- 
mer solution. Fig. [2]shows simulation results for 5R\\{t), 
measured from the projection on the initial longitudinal 
axis, and for fh{t) (proportional to the sum of the spring 
displacements from their equilibrium position). Param- 
eters were chosen such that the initial extension is close 
to full stretching (i?||(0) « 0.97i in all cases). While 
the universal scaling for longer times is evident (the ap- 
parent systematic offset in the "field" case arises simply 
because there is only one free end), substantial differ- 
ences between the scenarios for shorter times are clearly 
observable as well. 



3 



II. QUALITATIVE THEORETICAL RESULTS 

From Fig. [2] it is obvious that the time evolution 
of and fh{t) docs not obey simple power law 

scaling. S imp lifying approaches based on scaling argu- 
ments , clastic dumbbell models [H, [2^ , or quasi- 
equilibrium approximations [3, [131 have sometimes been 
used successfully for specific situations and parameter 
ran ges. In contrast, we employ a systematic formal- 
ism [35| based on the wormlikc chain model, which allows 
to generally account for the complex dynamics result- 
ing from different environmental perturbations. Here, 
we first present qualitative results for all four scenarios 
in order to illustrate their differences, and discuss exact 
analytical and numerical results in the next section. 

In the wormlike chain model [s^, semiflexible poly- 
mers are represented as inextensible smooth spacecurves 
r(s, t) of length L. Bending energy is proportional to 
the squared local curvature (9^r)^, such that in equilib- 
rium, tangent orientations are correlated over the persis- 
tence length £p = k/U-bT, where k is the bending rigidity. 
The initially straight polymer is supposed to be equili- 
brated at times t < Q and released at t = 0. After that, 
the longitudinal contraction is driven energetically up- 
hill via the creation of contour undulations by entropic 
forces. Considering the conservation of contour length 
due to the (near) inextensibility of the backbone bonds, 
these transverse wrinkles are conveniently referred to in 
terms of their excess contour length, or stored lengthy with 
an associated line density g. Mathematically, the inex- 
tensibility is enforced by the backbone tension f which 
counteracts stretching, and the creation of stored length 
is accompanied b y t he relaxation of tension. The the- 
ory of Refs. [3^, l39l | relates the tension f{s, t) to the 
stored length density q{s, t), based on the weakly-bending 
limit of small contour deviations from a straight line. 
In practice, this can easily be realized by choosing the 
control parameters /pro, E or v, or 7 sufficiently strong 
(as in typical experiments, see Table Hill below) . or the 
quenching factor 9 sufficiently large. It also justifies the 
free-draining approximation, where hydrodynamic effects 
are captured by anisotropic local friction coefficients Cj^ || 
(per length) for transverse/longitudinal friction, respec- 
tively [40[. However, ordinary perturbation theory is ap- 
plicable only for late times < ^ i^, ~ (\\L^ /(kBTi^) jSSj . 
because to lowest order it allows only a linear spatial de- 
pendence of g and / and neglects longitudinal friction 
forces [S^l- Further, except for quite stiff filaments with 
£p > L, the time is usually larger than the filament's 
longest relaxation time ~ C||£pL^/fcBr [4l|. which 
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within our approximations is given by the Rouse time 
of a polymer with Kuhn length 2£p. Nevertheless, with 
an improved formalism ^35i, i39] including nontrivial spa- 
tial variations in / and g, the conformational relaxation 
at times t ^ <h can be analyzed even for quite flexi- 
ble polymers. This leads to the remarkable insight that 
weakly-bending polymers constitute self-averaging sys- 
tems: the small stochastic fluctuations average out along 
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FIG. 3: (Color online) Schematic representation of conforma- 
tional (left) and tension relaxation (right) in a "force" setup. 
For t = 0, the filament is equilibrated under the force /pre- 
As the ends are released, the contour coils up in two growing 
boundary layers of size £|| (t) where the tension relaxes. At 
t = t^l, the dynamics crosses over from the propagation to 
the relaxation regime and the tension relaxes to zero. 



the contour and the coarse-grained tension dynamics fol- 
lows from the deterministic relation 



d!f = -C\\dt{g) 



(4) 



where the overbar denotes a (local) spatial average that 
produces effectively an ensemble average (denoted by 
( • )) [l^- Driven by tension gradients, stored length 
propagates subdiffusively from the filament's ends into 
the bulk — limited to boundary layers of size (t) by lon- 
gitudinal solvent friction. In more intuitive terms, the 
filament starts to "coil up" first at the boundaries, and 
only later in the bulk, see also Fig. [31 

In general, (g) is a nonlinear functional of /, see 
Eq. below for a detailed expression. Exact analytical 
results for the boundary layer size (t), the bulk tension 
fh{t), and the change in projected length <5i?||(t) will be 
obtained as leading-order results of a systematic asymp- 
totic expansion of Eq. Q in the next section. However, 
the scaling of the dominant part dR\\(t) of SR\\{t), which 
is independent of boundary conditions and an effectively 
deterministic quantity, can be found from a simple di- 
mensional argument: the change in end-to-end distance 
equals the amount of stored length p^y that has been cre- 
ated in the boundary layer. On the scaling level, Eq. ^ 
reads fh/£f\ — Ciip/^j s-^d '^g obtain 



Cii^l 



(5) 



Note that the change SRoit) in the gyration tensor's 
largest ei genv alue, which is frequently identified with 



SR\\ [211, 1421, obeys a different scaling law SRq — 

tfh/{C\\L) for short times t < t| Q. Additional sub- 
dominant contributions to (5i?|| (t) from end fiuctuations 
will be analyzed in section Hill 

Fig. m summarizes the scaling results of section HIJ for 
^11 (i), fhit), and dR\\{t) in various intermediate asymp- 
totic regimes, which are separated by different crossover 
times that have been matched by an appropriate choice 
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FIG. 4: (color online) Qualitative results. Asymptotic scaling laws for the boundary layer size lp{t) (a), the mean bulk tension 
,fh{t) or the bulk stress (b), and for the change in projected length SR\^{t) (c) in "force" (solid), "field" (dot-dashed), "shear" 
(dot-dot-dashed), or "quench" (dashed) setups, respectively. 



of the respective control parameters for better compar- 
ison. Clearly, the time t| = C\\L^[kYiT / {i^,fl^J\'^l^ is of 
key importance since it separates scenario-specific and 
universal relaxation. To understand the origin of the 
differences for times t <C we will consider the differ- 
ent scenarios separately before we address the universal 
regime t\. 

"Force" setup. After the stretching force has been 
shut off, the polymer starts to build up contour undu- 
lations driven by thermal noise. These transverse undu- 
lations appear first in growing boundary layers of size 
£|| [t) <IC L near the ends: assuming an inextensible back- 
bone, the immediate creation of undulations in the bulk 
would require the ends to be pulled inwards against longi- 
tudinal solvent friction with a force exceeding the actual 
backbone tension. This phenomenon of tension propa- 
gation ends after a time defined via -^[[(tf^) = L, 
where the boundary layers extend over the whole polymer 
length, see Fig. [3l A more detailed analysis ^] shows 
that the longitudinal relaxation depends on whether the 
internal tension (initially equal to the stretching force 
/pre) represents a relevant perturbation to the transverse 
conformational dynamics. The latter undergoes a dy- 
namic crossover from a bending-dominated regime with 
£|| cx t^/^ 01 for the shortest times t ^ tf to a tension- 
driven regime with oc t^/"^ [3l for longer times t ^ tf , 
where the dynamics becomes inherently nonlinear. Here, 

U = C-L^'B^^^p/Zpre is & crossovcr time that obeys ti <^ i| 
for reasonably large prestress /pre (see Table IIIII below) . 

"Field" and "shear" setup. Similarly, one can define 
a force equivalent in "field" or "shear" experiments and 
corresponding expressions for if and the hydrody- 
namic equivalent of /pre is simply the total Stokes friction 
C\\vL in a homogeneous fiow, and C||7-^^ is the longitu- 
dinal friction in an extensional shear flow. Flow condi- 
tions may straightforwardly be recast into the equiva- 
lent language of external (e.g., electrical) fields. How- 
ever, complicated counterion effects [45|, \^ pre- 
vent the quantitative prediction of the equivalent elec- 



trophoretic field strength E for typical experimental re- 
alizations. Although field-type perturbations induce dy- 
namic crossovers at if similar to the "force" case, the 
change SR\i {t) in projected length increases always lin- 
early with time. This can be understood by a simple 
change in perspective: the polymer's ends are pulled in- 
wards by an approximately constant bulk tension, i.e., 
with roughly constant velocity. This corresponds in 
the frame of reference of the ends to an external flow 
field. The resulting friction forces are properly balanced 
and the initial polymer conformations are already equi- 
librated under such a flow field in "field" and "shear" 
setups, in contrast to the "force" scenario. Tension prop- 
agation is therefore not a dominant effect in the former 
(Eq. ([5]) applies with £|| = L because of the large-scale 
spatial variation of the tension), and the constant drag 
gives (5^11 (X t. 

"Quench" setup. Here, finally, there is no external 
force scale and therefore no dynamic crossover (^y cx 
^1/8 fQj- I ^ ^1)^ although the parameter combination 
kBTipO'^ / plays the role of /pre in the crossover time 

. The tension in a quenched filament is produced solely 
by the suddenly increased thermal noise from the envi- 
ronment (if the external temperature increases), or by 
the suddenly higher "sensitivity" to this noise (if the 
quenching is achieved by a sudden drop in bending rigid- 
ity). Hence, its magnitude depends on the "mismatch" 
between the current conformation and an equilibrium 
conformation corresponding to the current environment. 
Therefore, the quenched filament can relax tension even 
in the bulk, by reshuffling stored length between long and 
short wavelength modes in a way similar to the mechan- 
ical stress relaxation in buckled rods [i^ , while the bulk 

tension stays constant for i <C t| in the other setups. 

Universal regime. At long times i 3> t|, when the 
tension has propagated through the filament, the dynam- 
ics enters the universal regime of homogeneous tension 
relaxation [35| . Contrary to previous assumptions [20| . 
longitudinal friction may not generally be neglected, but 
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dominates the dynamics in this regime. The tension has 
a nontrivial spatial dependence, but it can for asymp- 
totically large forces be treated as quasi-statically equi- 
librated [13, UHl- The characteristic universality of the 
long-time relaxation is then simply a consequence of the 
right hand side of Eq. ([4]) being independent of initial 
conditions: g ~ [ksT/ (4£p/)]^/^, where we have used the 
(static) force-extension relation for wormlike chains [i^ . 
This asymptote readily implies by Eq. Q the scaling 
/ oc and by Eq. ([5]) the characteristic i^^'^-growth 

of SR»(t}, which has indeed been observed in experi- 



ments [22[. As an aside, we note that < ipi for stiff 
polymers with L < £p; the adjoining regime of alge- 
braic relaxation for times t^, <C t shows a i^/^-scaling 
in SR^it) P,[4i]. 

Let us finally comment on the joint limiting scenario: 
the exactly straight initial conformation (as in Ref. [IH). 
Not only it is quite artificial from a theoretical and ex- 
perimental perspective, it also appears to be ambiguous, 
since we could let /pio 00 in one of the scenarios in- 
volving external forces as well as — > 00 for the "quench" 
case. Although i| — > in both cases, so that only the 
universal regime survives and the ambiguity is limited 
to t = 0, it gives rise to observable effects as soon as 
one takes into account some microstructure corrections 
important for real experimental systems and simulation 
models [soj . 

After this qualitative discussion of the relaxation dy- 
namics, we will now present a systematic derivation and 
analysis of Eq. ([4]), resulting in exact growth laws in 
various intermediate asymptotic regimes for the observ- 
ables introduced above. While Ref. only covered the 
"force" case, we now obtain results for the other scenar- 
ios as well, and include a quantitative analysis of different 
boundary conditions. 



III. QUANTITATIVE THEORETICAL RESULTS 

Starting point for our calculations is the wormlike- 
chain Hamiltonian 



(6) 



where the backbone tension f{s,t), a Lagrange mul- 
tiplier function [5l|, takes care of the local inextensi- 
bility constraint r'(s)^ = 1. The equations of mo- 
tion for the contour result from balancing clastic forces 
—dH/dr with stochastic noise ^ and anisotropic vis- 
cous local friction forces Ci^tf — u] with friction ma- 
trix C = [Ci + C||(l ~ r'r')] and a velocity field u 
of the solvent. The friction coefficients (per length) are 
(^1, = with C ~ 1/2 and C_L ~ 47rr;/ln(L/a) 0, 
where a is the backbone thickness. In all of this sec- 
tion, we set = K = I for simplicity. This makes 
time a length^ and the tension a length"^. Note that 
K = k^Tip is kept constant in the "quench" scenario. 



Our approach exploits the weakly-bending limit. Pa- 
rameterizing the contour r = {r±,s — r^^)'^ in terms of 
small transverse and longitudinal displacements from the 
straight ground state, this means that r'l = 0{e) 1, 

with s = /i^y^/^p for "force" setups (and /pre replaced 
by its equivalents in "field" or "shear" scenarios) and 
e = L/{6£p) for "quench" setups, respectively. Up to or- 
der e, the equations of motion for the contour in absence 
of external forces and for m = read: 

ftri = -rr + (M)' + a (7a) 
Ca^rii + (1 - Qr'^dtr^ = -r[[" - /' + (/r[| )' + (7b) 

Because in the weakly-bending limit the transverse con- 
tour fluctuations arc correlated on much shorter length 
scales than the longitudinal (= tension) dynamics, wc 
can formally introduce "fast" and "slow" arclength co- 
ordinates for the small-scale transverse and large-scale 
longitudinal dynamics, respectively [s^. Taking a local 
(with respect to ^y) spatial average over the small-scale 
fluctuations (denoted by an overbar) leads to closed equa- 
tions: 



dU=~C9t{g). 



(8a) 
(8b) 



■^r'^ is the 



The longitudinal part Eq. (|8bl) . where g = ^r'l 
stored length density, follows from the self-averaging 
property of weakly-bending polymers: the spatial coarse- 
graining effectively generates an ensemble average. The 
transverse part Eq. ([8a|) contains a locally constant ten- 
sion / (its slow arclength dependence obtained through 
Eq. (j8bp is adiabatically inherited), and can be solved in 
terms of appropriate eigenmodes w<j(s) with eigenvalue 
—q^{q^ + f{t)) via the response function 



(9) 



Using the noise correlation {$±{k,t)^±{q,t')) = 
A£p^Sk.qS{t — t'), we evaluate the expectation value 
(ir^). The different preparation mechanisms discussed 
in the main text constrain the polymer only for t < 0. 
Including the initial conditions f{s,t < 0) = /o(s) and 
£p{t > 0) = Ip ^ £p{t < 0)/9 gives for the stored length 
density 



(9) 



■E 



q'0{q^ + Ms)) 



+ 2 / dt' xliq;t,t') 



w'hs) (10) 



As only the spatially averaged stored length density {g) 
enters Eq. (j8bp . we decompose w'^{s) into a spatially con- 
stant and a fluctuating part Cq{s) (the latter will average 
out upon coarse-graining): 



<is) 



L 



[l + c,(s)] 



(11) 
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Taking the continuum limit L — > cxj and integrating 
Eq. (|8b[) over time, we find that the integrated ten- 
sion F{s,t) — /gdi'/(s,t') obeys the partial integro- 
differential equation 



dq_ 

TT^p 



i-xi(g;^,o) 

0{q' + /o(s)) 



2q' I dt'xl{q;t,t') 
10 



(12) 



From solutions to this equation in different intermediate 
asymptotic regimes presented in the next subsection, we 
will then infer growth laws for the two observables. 



A. Asymptotic results for the tension 

1. "Force" setup 

This scenario with 6 ~ I and the initial and boundary 
conditions 



/(s,f < 0) = /pi-c and 
/(0,t>0) = 0, /(i,<>0) = 0, 



(13) 



is identical to the "release" -scenario which was thor- 
oughly analysed in Ref. [4l|. We will briefly sketch this 
analysis in order to motivate its application to the other 
setups. From the response function Eq. we get the 
asymptotic scaling for the wave number Q of the mode 
that relaxes at time t: 




ifFVt<l ("linear") 
ifFVi>l ("nonlinear") 



(14) 



Examining Eq. ([5a|) . one infers that in the first case the 
tension contribution Q^F is small compared to the bend- 
ing contribution Q'^t and can be treated as perturbation 
on the linear level. Since the magnitude of the tension 
is determined by the prestretching force, F ~ /pre^, this 
asymptote, called "linear regime" , can also be formulated 
as t <^ tf with tf = In the second case t :s> tf, the 

bending contributions are subdominant which leads to 
different "nonlinear regimes" . 

Linear propagation (t <C tf). We perform an expan- 
sion [H^l of the right hand side of Eq. ([72]) with respect 
to the integrated tension F and to the force /pro: 



+ 2F{s,t)-Aq^ / dt'J^(s,i')c"^'''^*"*'^ 

"'0 



(15) 



Using the Laplace transform F{s,z) = C{F{s,t)}, this 
reads 



dq 



2/, 



pre 



z{z + 2q*) 



+ F{s,z)-^ 



2z 



(16) 

which, after performing the q-integral, reduces to: 



A2 a^F = F 



fpr 



pre 



(17) 



Here, A(z) = 2^/*^(£p/C)^/^z"^/^ is a dynamic length 
scale denoting the size of spatial variations in F(s,z). 
If L ^ A, the solution to Eq. p7|) varies only close to 
the boundaries, as it is characteristic for the propagation 
regime. Near s = (and correspondingly near s = L), it 
simplifies to 



F{s,z) « 



l-e 



-s/A 



which can be backtransformed |4li to 



i^(s,i)=/prct[l"0(s/^|[(t))] 



(18) 



(19) 



where 0(^) k. exp[— 2~^/^^/r (15/8)] is a scaling func- 
tion that depends only on the ratio ^ = s/f||(t). The 
length scale A is directly related to the boundary layer 
size l\\{t) = (l^llfl'^t^l^ P, El, m, and the require- 
ment L 3> A translates into t <C i| . 

Nonlinear propagation (tf <^ t <ti In the nonlin- 

ear regime, the QH bending contributions are small com- 
pared to the tension terms Q^F ii F'^/t ^ 1. This results 
in x±{Q': t') being finite only near t' w t, see Eq. ((9)). We 
can therefore hnearize f'(s,t)—F(s, f') « [dtF{s,t)] {t—t') 
in the exponent. The i'-integral in the second term of 
Eq. is readily performed ^41j : 



^dq_ 



C 3- 



dg 



+ /pre 

1 



dtF 



q" + /pr. 



dtF 



(20) 



In the second line, we let x± ^ because F'^/t ^ 1. 
This indicates the underlying "quasi-static" approxima- 
tion: the relevant modes have already decayed and the 
tension is quasi-statically equilibrated. Taking a time 
derivative gives [13, EH 



4€p/3/2- 



(21) 



Inspired by the result Eq. (|T9|) , we expect a scaling form 
f{s,t) = fpio'piO with ^ ~ s/£||(t) for the tension. In- 
serting it intoEq. (ini) gives £||(i) = (fp/O^/Vp™ i^/^ 0, 
m and 



(22) 
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with the boundary conditions ip{0) = and d^(p{^ 
oo) = 0, i.e., we neglect the presence of the second 
end, where the situation is correspondingly, and assume 
just a flat profile in the bulk. Numerical solutions to 
this equation have been shown in Refs. EH and give 
(p{£_ oo) = 1 as expected and ip'{0) ~ 0.62. The prop- 
agation regime ends at t\ ~ fp^V^ / ip. 

Homogeneous relaxation (t\ <C i <C t^)- After the 
tension has propagated through the filament, it is no 
longer constant but expected to decay. But as long as 
F"^ /t ^ 1 still holds, we can use Eq. (PT|) . Hence, we try 
the separation ansatz f{s,t) = g(t)h{^) with ^ = s/L, 
which gives [4l| 



git) 



2/3 



and 



h" 
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-1/2 with /i(0) h{l) ^ 0. 



(23) 



(24) 



The almost parabolic profile h{^) is characterized by [41| 



h'{0) = 12-l/^ Ml/2) = 



I 2/3 



(25) 



Using Eq. (gS]), wc find that the condition F"^ jf^ 1 is 
violated for t >ti, = /£p, which is already larger than 
if £p < L. Hence, this regime lasts until the weakly- 
bending approximation breaks down near the ends due 
to the onset of the stretch-coil transition. 



2. "Field" setup 

For hydrodynamic and/or electrophoretic forces, we 
find from the longitudinal equation of motion Eq. (j7bp a 
corresponding non-uniform initial tension profile /(s, t < 
0) = g{L — s) with g ^ (v for flows or g (x E for an 
electric field, where the generally unknown prefactor is 
some combination of electrophoretic and hydrodynamic 
mobility. This linearly decreasing prcstress would in prin- 
ciple lead to an additional term f'r'j_ in Eq. (|5a|) . and the 
corresponding eigenfunctions would be very complicated. 
However, because large scale tension variations are irrel- 
evant for the short wavelength transverse dynamics, we 
can ignore this term by consistently exploiting the scale 
separation which allowed the derivation of Eq. ([8]), and 
use Eq. (|12p with the initial linear profile fo{s) ~ g{L—s). 
The polymer is supposed to be grafted at s = and to 
have a free end at s = i, i.e., the boundary conditions 
are 



/'(0,t>0) = and 



/(L,i>0) = 0. (26) 

Identifying the force equivalent f*^.^ — gL, we expect a 
linear regime for t ^ t[ with t[ = fp~^ and a nonlinear 
regime for t ^ ti, governed by the respective asymptotic 
differential equations from the "force" case. 



Linear propagation (t <C tf). Linearizing Eq. (jl2p in 
F and /o(s) and performing a Laplace transform as in 
Eqs. p^^TT)) . we arrive at the solution 



F{s, t) = gt{L -s)- (i)t0(s/^|| (t)). 



(27) 



,3/8 



oxp[ 



£(17/8) 



With ^ii(t) = {epKY'H^'^ and m . 

We find a boundary layer at the fixed end where the ten- 
sion relaxes from its initial value gL only by the small 
amount . Near the free end, at s = L, we have 
F{L,t) = and F'{L,t) = —gt without any algebraic 
correction terms. Hence, because the tension at the free 
end is already very small and the contour does not fur- 
ther coil up, there are no boundary layer effects which 
would give relevant deviations from the linear drift to- 
wards t he g rafted end, in contrast to what has been found 
in Ref. [l7| . 

Nonlinear propagation (t^ t <^ t\). The assump- 
tion F"^ /t ^ 1 leads again to Eq. (pT)) except for very 
small regions near the free end where /o(s) in the denom- 
inator of the first term of Eq. (|12p is almost zero. Corre- 
sponding to the linear case Eq. ((27)) . we assume that the 
tension deviates only near the fixed end from its initial 
value /o(s). Hence, we insert f{s,t) — /o(s) — g£\\{t)'p{0 
with ^ = s/e\\{t) and £||(t) = {ep/CY^^igLf/H^/^ into 
Eq. (pij). and expand about /o(s <C L): 



1 



(28) 



with (p'{0) = —1 and (p{£, oo) = 0. The solution can 
be given in terms of the complementary error function: 



^(e) = ^c-^'/i6-ecrfc(e/4). 



(29) 



The propagation regime ends at rj^ = L" /{£p(gL)^'-). 

Homogeneous relaxation (tf^ "C i <C tR). The separa- 
tion ansatz f{s,t) = g{t)h{S,) with ^ = s/L in Eq. ((2T|) 
gives g(t) = {C^L"^ / [Ipt))"^/^ as in Eq. ^ and the spatial 
function h solves 

h" = -Ih-^'"^ with h'{Q) ^ 0, h(l) = 0. (30) 

We find the following characteristics 



M0) = (A)2/3^ /.'(I) = 6-1/3 

The condition F^/t:$> 1 holds until t = t^ . 



(31) 



3. "Shear" setup 

In this scenario, the equations of motion Ci^t''' — u] = 
—5Ti/6r+^ are modified in the presence of an extcnsional 
shear flow field u = j{~r±, s~r^^ — L/2)-^, where 7 is the 
shear rate. To lowest order in e, and in the stationary 
state, we obtain from Eq. ([Tb]) 



(32) 
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As before, we treat this non-uniform tension profile only 
as large-scale variation and use Eq. with the initial 
and boundary conditions 



(33) 



/(s, t < 0) ^Cis {L - s) and 
f{0,t>0)^ f{L,t> 0)^0. 

As in the "field" case, the time tf = /p^^ with the 

force equivalent f*^^ ~ CiL"^ denotes the linear-nonlinear 
crossover. 

Linear propagation (t <^ tf). Here the solution to 
Eq. ini) reads 

Fis, t) = iC7is(L -s)- -YJ^ [1 - ^('/h W)] , 

(34) 

with = (^p/O^/^ti/® as before and 0(0 ~ 

exp[-2-3/8r (9/4) C/r (17/8)]. We find two smaU bound- 
ary layers at the ends where the tension is slightly smaller 
than initially. 

Nonlinear regime (t{ <^ t <$:i The nonlinear 

regime tf ^ t <C t| for the "shear" setup is quite pe- 
culiar: if we try (similar to the "force" and "field" case) 
a scaling ansatz /(s, i) — ^C,'^[s{L — s) + L£^^{t)ip{s/£^^{t))] 
or similarly, we get i\\{t) ^ t^. This unusual result could 
be explained by the fact that the "prestress" f{s,t < 
0) w ^C-ysL, which is responsible for the scaling of £|| in 
this regime [52| , grows linearly with the distance from the 
ends. However, we dot not get any physically meaning- 
ful differential equation for ip under the boundary condi- 
tions Eq. ([33]) . We conclude that there is no propagation 
and no observable boundary layers. Looking for a solu- 
tion spanning the whole arclength interval from to L 
instead, we insert into Eq. (|2ip an expansion of the form 



f{s,t) = ^CiL' 



MO + ^MO + o{{t/tl) 



(35) 



with ^ = s/L and t[ = (L^ /[ip{('^L^f^^]. Solving the 
resulting differential equations for successive powers of 
(t/t^l) gives the leading order terms 



3/2 



(36) 



In contrast to the propagation forms f{s,t) ~ ip{s/£^^{t)) 
of the other scenarios, we now get self-similar and spa- 
tially invariant tension profiles. This can probably at- 
tributed to this specific initial condition which allows 
for self-similar relaxation. Further, to linear order in 
(t/t^l) we do not obtain algebraic corrections to the linear 
growth law of SR\\ {t), because d^(pi{0) = 0. Higher-order 
terms in the expansion (as far as they are analytically 
tractable) turn out to be ill-behaved near the ends. 

Homogeneous relaxation (i(:| <C t <C t^,). The subse- 
quent regime of homogeneous tension relaxation is ex- 
actly equivalent to the one of the "force" case (see the 
respective boundary conditions Eqs. (fT3l [33]) ). and the 
results Eqs. (^51 1^5)) apply here as well. 



4. "Quench" setup 



This scenario, with the initial and boundary conditions 



f{s,t <0) = 
£p{t <0) = ei 



and 



/(0,t>0) = /(i,t>0) = 
ipit > 0) = £p 

has been introduced as "£p-quench" in Ref. [35[. In con- 
trast to the scenarios discussed above, we lack a quan- 
tity providing a force scale, and the tension attains a 
simple scaling form in the propagation and relaxation 
regimes, i.e., there is no linear-nonlinear crossover. How- 
ever, this scaling form still strongly depends on the value 
of 6. Physically the limit 9^0 corresponds to sud- 
denly switching off thermal forces for a thermally equi- 
librated filament, hence purely deterministic relaxation 
(see Ref. [43^]). A small quench « 1 will not in- 
duce strong tension, but the limit 9 00 describes 
the scenario of a completely straight contour that equi- 
librates purely under the action of stochastic forces, 
and the resulting tension may be very large at short 
times. In this case, similar approximations as employed 
when discussing the nonlinear regime in Sec. IIII A II can 
be justified. Assuming F'^/t ^ 1, the response func- 
tion X-L('?;tjt') from Eq. ([9]) is finite only near t' w 
t which suggests the linearization F{s,t) — F{s,t') w 
[dtF{s,t)]{t — t') in the exponent. Performing the t'- 
integral in the second term of Eq. ((T2)) yields 



'dq_ 



i~xliq;t,o) i-xi(<z;t,o) 



dtF 



In contrast to Eq. ([^0)1 . we may not set x± ~^ in the 
first term, because this would produce an IR-divergence. 
But we can neglect the bending contribution qH in the 
exponent of the first and set x± ~* only in the second 
term. This gives 



'dq_ 



1 - e-^^"^ 



1 



9q^ 



dtF 



IF ^ {dtF) 



-1/2 



(38) 



Propagation (t <C t\). Inserting the scaling ansatz 
F{s,t) = 9t^/^^{s/£\\it)) with (t) = (£p/C)_i/203/4ti/8 
removes the parameter dependence in Eq. 



-1/2 



(39) 



Boundary conditions arc '/>(0) = and d(_4>{£^ 00) = 
0. From a numerical solution we obtain ^(^ 00) ~ 
0F/2 as expected from Eq. ^ and 0'(O) « 1.44. The 
assumption F"^ /t ~ 6* ^ 1 is justified for all times t ^ 

t\ = Lyiip% 

Homogeneous relaxation ft}]^ ^ t ^ tj^J. Because 
we expect universal long time relaxation in the strong 



quenching limit 6 ^ 
the ansatz F(s, t) 
Eq. 



• oo similar to the "force" -case, we try 



(40) 



Because now f{s,t) = dtF{s,t) = g{t)h{£,) with = 

{(L^/iptf/^ and = (/-(O/S as in Eqs. (QllMl), this 
regime of homogeneous relaxation is identical to the one 
in the "force" -case. The condition /t ^ 1 holds until 




than if ^ i. 



which is usually already larger 



B. Results for pertinent observables 

The maximum bulk tension fh{t) = f{L/2,t) (in the 
"field" case, we prefer to use the grafting force fg{t) = 
f{0,t)) can be obtained directly from the tension pro- 
files computed in the preceding section. The change in 
end-to-end distance (t) follows from a simple formula: 
With the sign convention used before, this change has to 
equal the total amount of stored length that has been cre- 
ated. Hence, we integrate the ensemble averaged change 
in stored length density dt {g) {s, t) over s and t: 



8R\\{t) 



d.S / dt'dt{Q){s,t') 

Jo 



(41) 



Defining {g) = (g) + {g°) in Eq. pTT)) from the decompo- 
sition Eq. dll]), we obtain (5i?||(t) ^ SR\\{t)+S^{t). The 
first part obeys the deterministic Eq. (|8bp : 



Jo Jo [^^) 
= ~C' [dsFiL.t) ~ dsFiO,t)] . 

It accounts only for the "slow" coarse-grained tension 
dynamics but neglects subdominant and stochastic con- 
tributions (5-R|j(t) from "fast" fluctuating boundary seg- 
ments analyzed in the next section. Results for SR\\ (t) 
and fh{t){fg{t)) are summarized in Tables HI and [U 





t < t{ 


tc <t^t{ 


4 <^ 


"force" 




fl/4 1/2 


( ISLt A 


r(i5/8)ci/24/2 




\ p / 


"field" 


gt/C 










/ \ 1/3 


"shear" 




■ytL 


/ ISLt \ 


"quench" 


2.88 


I ISLt \ 

[ c^i ) 



TABLE I: Asymptotic scaling laws for the change in projected 
length (5_R|| (t) from Eq. (|42p for the different setups. Units 
have been chosen such that k, = (± = 1. 





t < tc 






"force" 


/pro 


( 3CL2 yr^ 

\ 128£pt J 


"field" 


gL[l-Om/L)] 


\ 32£pt J 


"shear" 




^(jL^[l - 0{t/t[)] 


( 3CL^ 
\ 128fpt J 


"quench" 




( 3CL^ yrs 

\ 128£pf J 



TABLE IL Asymptotic scaling laws for the maximum bulk 
tension /b(t) = f{L/2,t) (grafting force fg{t) = /(0,t)) for 
different setups. Units as in Table IT] 



conditions for the contour r(s) is only subdominant but 
still non-negligible. While the "bulk contribution" SR^^ 
is independent of boundary effects and dynamically self- 
averaging, this stochastic dependence is accounted for 
by an additional "end contribution" 6Rp which stems 
from the oscillating term Cg(s) in Eq. and decays 

rapidly on much smaller length scales than that of ten- 
sion variations [HI. We may therefore evaluate this 
part at the boundaries under zero tension (i.e., using 
X±iq;t,t') = c-^'^f*-*') instead of Eq. 



ds 



dq q^-e{q^+Us)) 



x(l- 



/o(^)) 



(43) 



C. Boundary effects 

Because our theory apphes to times t <^t^ long before 
the relaxation of long-wavelength modes becomes rele- 
vant, because Eq. ^ results from a coarse-grained de- 
scription that averages over small-scale fluctuations, and 
finally because projecting the end-to-end distance onto 
the longitudinal axis suppresses some end effects (4l| . 
the dependence of 6R\\ = 5R\\ + (5i?|j on the boundary 



Consistent with this simplification, we approximate the 
Wq{s) by eigenfunctions of the biharmonic operator 
(see Ref. 53|). Again exploiting the scale separation in 



this integral over the rapidly fluctuating term Cg(s), we 
use only the spatial average /b of the slowly varying pre- 
stress /o(s). This means /pre — »■ /b = j2^^^'^ 
"shear" case, and /pio ^ /b = ^gL for the "field" setup. 
Because there is only one free end in the latter, the con- 
tribution to 6Rn is one half of the "force" result. 
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Free ends. If w'^ = w'^ 



at s = 0, L, we use 
sin qL + sinh qL 



- {sin qs + sinhgs) 



^/ZlcosqL — coshqL 

+ cos qs + cosh 



(44) 



where g is a solution of cos cosh = 1. For t t^i, 
the q-integral in Eq. is dominated by short wave- 
length contributions, and for the relevant asymptotically 
large modes the s-integral over Cq{s) reads 



/ dscg(s) 
Jo 

Hinged ends. For Wq = 

Wq{s) : 



6 



(45) 



at s = 0, L, we obtain 



J- sm qs 



(46) 



with singL = 0. In this case, dscg(s) = 0. 

Clamped ends. Here, Wq = w'^ = ai s = 0, L, and 
we have 



1 



sin qs — sinh qs 

cos qL — cosh qL 
sin qL + sinh 



(cos qs — cosh qs) , (47) 



with cos qL cosh qL = 1 and dscq{s) ^ — 2/q. Up to a 
prefactor, the contributions for clamped ends is identical 
to the one for free ends. 

Torqued ends. If w'q — w'q = at s = 0, L, the Eigen- 
modcs are 



(48) 



Wg(s) = \ J- cos qs 



with sin qL = and dsCfy(s) = 0. 

Using the asymptotic limit of Cq(s)ds, we evaluate 
the end contributions in the free ends situation of our 
setups: 

"Force" setup. Here, we obtain 



(49) 



where 
G(0) 



dfc- 



1 



'o fc3(fc2+0) 

^[c-2^^(7rcrfiy2^~Ei(202)) 



(50) 



with erfi the imaginary error function, Ei the exponential 
integral and 7e ~ 0.577 Eulcr's constant. The asymp- 
totic behavior is 



3/prc^ ,„ o J-2 



5Rm 



TT^p 



27r€„ 



t < /pre (51a) 

t » f-l (51b) 




"Quench" setup. For free ends, we obtain from 
Eq. (021): 



(52) 



As anticipated [4l|], the end contribution (5i?jj is zero 
for hinged and torqued ends, while it leads to an addi- 
tional reduction of for free ends and has a lengthening 
effect for clamped ends. Note that indeed <C 5R\\ is 
only subdominant for t <C tR, but the quantitative rel- 
evance of this contribution becomes evident when it is 
directly compared to numerical solutions of Eq. p2)) and 
simulation data in non-asymptotic regimes. 



IV. COMPARISON TO SIMULATION DATA 

The asymptotic scaling laws of Fig. 3] are derived in 
the limit /pro ^ oo, and the difference to the numer- 
ical solutions gets smaller than 20% only for /pro ^ 
lO^^fceL^p/L*. While this can easily be realized in ex- 
periments, for instance on DNA (cf. Table IIIII below), 
it is not possible in simulations due to the usual trade- 
off between computational efficiency and accuracy. For a 
comparison between our theoretical results and the simu- 
lation data of Fig. [21 we therefore compute the bulk part 
5R\\{t) and fh{t) using numerical solutions to Eq. (fT^ 
as described previously [s^ . Fig. [5l shows simulation 
data and analytical results for all four scenarios, using 
/i = (C±&)^^ and = C|| to relate the (isotropic) mobil- 
ity of the beads in the simulation model to the anisotropic 
friction coefficients per length of the continuous wormlike 
chain used in our theory. The bulk contribution 5R\\ (t) 
(dashed lines), while having the correct qualitative be- 
havior, underestimates the contraction by as much as 
50%. Including end fluctuations with (5i?jj(i) (solid lines) 
gives results for 5R\\ (i) that are slightly overestimated for 
longer times. This could be caused by possibly oversim- 
plifying approximations made when evaluating Eq. (|43p . 
or by a gradual breakdown of the weakly-bending limit. 

In the "quench" case, we observe a strong deviation 
between simulation and theory for short times, both in 
5R\\ and /b. While /b cx: diverges as i — *■ in 

our theory, the actual tension in the simulations is fi- 
nite. This is due to the extensible backbone of a head- 
spring chain: the tension follows the sudden change in 
environmental conditions only with a temporal delay re- 
lated to the finite propagation speed of longitudinal back- 
bone strain. Because now /b is smaller than predicted, 
the contraction 5R\\[t) is also reduced, see the scaling 
relation Eq. It is, however, possible to include a 

finite extensibility correction in Eq. Because this 

nontrivial extension is only marginally relevant for the 
present discussion, which is focused on differences in the 
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FIG. 5: (Color online) Comparison between theory and simulation data for a "force" (a), "field" (b), "shear" (c), and "quench" 
(d) scenario, respectively, for (5_R||(t) and fh{t) (insets). Note that the bulk contribution SR^^{t) from Eq. (|42|) (dashed lines) 
underestimates the contraction, which is corrected for by including the end contribution 5R°^{t) from Eq. (|43|l (solid lines). In 
the simulation of the "quench" scenario, the tension follows the sudden change in temperature only with a delay related to 
the longitudinal propagation of backbone strain, which leads via Eq. (O to a different scaling of at short times. This can 
be accounted for by a correction term (dotted line) including the finite backbone extensibility of the bead-spring model [50| . 
Simulation data as in Fig. (2] 



relaxation dynamics from an initially straight conforma- 
tion, we present a detailed discussion elsewhere [lOl , and 
merely show the corrected results for {t) and fh{t) for 
the "quench" case in Fig.O^d) (dotted lines). The anal- 
ysis of this correction term allowed to choose parameters 
such that our results are not affected by microscopic de- 
tails for the other setups, but this was not possible for 
the "quench" case with its singular short-time behavior. 



(a) 


DNA 


actin 


(b) 


DNA 


actin 


tf 


10" ' s 


10"" s 


/c 


0.08 pN 


0.2 fN 


til 

''L 


0.05 s 


0.003 s 


Vc 


27 A*ni/ s 


72 nm/ s 


ta 


^ 6s 


^ 10s 


7c 


6.6s-^ 


0.02 s"^ 



TABLE III: (a) Characteristic time scales (given a force of 
/pro ~ 2 pN) and (b) bounds on control parameters for typical 
DNA [3 (i ~ 20 ^tm, £p « 50 nm) and actin [sS] (i « 20 ^m, 
£p 17 fim) in solution with viscosity 77 ~ 10"'^ Pa s at room 
temperature. 



V. EXPERIMENTAL IMPLICATIONS 



Altogether, we now obtain good quantitative agree- 
ment between computer simulation and theory for all 
four setups and both observables over six decades in time 
without adjustable parameters. Having reliable theoret- 
ical control over the relaxation dynamics, we will now 
present quantitative estimates for the feasible choice of 
control parameters in experiments and a qualitative dis- 
cussion of the influence of some additional important ef- 
fects. 



A. Time and force scales 

In Table IIIII we have compiled numerical examples for 
the various time and force scales introduced above based 
on literature values for DNA [2^ and actin [s^]. In or- 
der to obtain sufficiently straight initial conformations, a 
conservative estimate for the control parameters /pro, 
and 7 requires them to be chosen by a factor 25 larger 
than the respective values /c, Vc, and jc- The quenching 
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strength 9 should be significantly larger than 9c = L/£p. 

The crossover times t| (x L'^/ fp(c and t{ cx /p^^ depend 
strongly on the adjustable quantities L and /pro; hence 
the time window of interest can be varied considerably 
between scenario-specific and universal relaxation. The 
algebraic relaxation ends at times near tR, for which we 
can give only a rough estimate as the unknown numer- 
ical prefactor is influenced by boundary conditions and 
hydrodynamic interactions and may substantially differ 
from unity (30l | . 



B. Onset of the stretch— coil transition. 

Since experiments are often performed using quite flex- 
ible polymers like DNA with £p <C L, the weakly-bending 
approximation will finally become invalid in regions near 
the ends, where the contour starts to (literally) coil up as 
the tension relaxes. Borrowing ideas from flexible poly- 
mer theory allows to derive scaling laws accounting for 
the onset of the stretch-coil transition. The stem-flower 
picture of Brochard-Wyart [l3| describes transient relax- 
ation processes of flexible polymers with Kuhn length a 
and friction coefficient (.. Entropic forces on the order 
of k^T / a arising in the bulk pull the ends inwards. Bal- 
ancing these forces with the associated friction gives the 
well-known scaling ~ [fceTt/i^a]-^/^ for the growth of 
"flowers" leading to an additional longitudinal contrac- 
tion. 

In the case of strongly stretched semiflexible polymers, 
this correction is negligible on time scales t ^ t^,. Here, 
the Kuhn segments arc of size €p, and their Rouse-like re- 
laxation after internal bending modes have equilibrated 
would generate flowers of size ~ L{t/t^y^^. How- 
ever, in the relevant universal regime of homogeneous 
tension relaxation (tl^ ^ t <C t^), the bulk tension 
/b ~ kBT(tB./t)^^^/£p (see Table HH is much larger than 
k-^T / Ip, and the ends are pulled inwards so fast that 
a flower of the above size would be too large for the 
resulting drag. Observing that the associated roughly 
parabolic tension profiles [4lj attain values of about 
k^T/ip within distances l'^ ~ L{t/tYif/^ from the ends, 
one easily confirms that this smaller value for the flower's 
size indeed restores the friction balance. Altogether, 
we find that for times t <^ the (t/tR)^/^-growtli of 
"flower" -like end-regions, where the weakly-bending ap- 
proximation breaks down, is subdominant against the 
(i/ipt) ^/^-contraction of the remaining weakly-bending 
part of the filament. Only at t ~ i^, our assumptions 
finally cease to hold and more appropriate models, for 
conformational relaxation as well as hydrodynamic inter- 
actions, need to be employed (see, e.g., Ref. [l8|)- In our 
simulation, we do not expect a pronounced stretch-coil 
transition because the number L/ £p is not large enough. 
We also checked that global rotational diffusion [i^l , ap- 
parently reducing the longitudinal projection of the end- 
to-end distance, can be neglected. 



C. Hydrodynamic interactions 

Finally, we want to briefly comment on hydrody- 
namic interactions. Their pronounced effects for strongly 
coiled polymers reduce to mere logarithmic corrections 
for relatively straight filaments [401. suggested pre- 
viously [1^, we speculate that these corrections can 
be summarily included via a phenomenological renor- 
malization L — > t^s in the friction coefficient C, oc 
ri/hi{L/b) [4^, where 77 is the solvent's viscosity and h 
the monomer size. Within our model, this has almost 
no further consequences than slightly shifting the time 
unit, cf. Eqs. @ and ([5]). An appropriate time rescaling 
compensates for changes in the friction coefficients, and 
could therefore easily be checked in experiments. The 
setups of Refs. [2^ [i^], where initially stretched DNA 
relaxes with one end attached to a wall and the other 
fixed to a bead, can easily be modeled within our the- 
ory by appropriately adjusting the boundary conditions 
for the tension. It turns out that spatial inhomogeneities 
of the tension and end fluctuations are suppressed and 
hydrodynamic interactions (primarily between bead and 
wall) are enhanced, such that a simple quasi-stationary 
approach describes the data very well. 



VI. CONCLUSION 

We have presented a comprehensive theoretical anal- 
ysis of the conformational relaxation dynamics of semi- 
flexible polymers from an initially straight conformation. 
Special emphasis has been put on a systematic investiga- 
tion of four fundamentally different realizations of "ini- 
tially straight". The sudden removal of the straighten- 
ing constraint leads in all cases to strong spatial inho- 
mogeneities of the filaments' backbone tension. Analyz- 
ing two exemplary and easily accessible observables, we 
found that for short times, when these nontrivial spa- 
tial variations are restricted to the boundaries, the re- 
laxation dynamics crucially depends on the actual ini- 
tial conditions: polymers pre-stretchcd with forces dis- 
play tension propagation effects, in contrast to chains 
straightened by fields or flows, and a quench leads to 
yet other effects. In the universal relaxation regime at 
longer times, the tension becomes quasi-statically equili- 
brated and independent of initial conditions, but its spa- 
tial inhomogeneity remains relevant. Additionally to the 
derivation of asymptotic growth laws, we extended the 
systematic theory of Refs. [s^ [s^ |4l| to include the sur- 
prisingly important influence of different boundary condi- 
tions. For non-asymptotic parameter values, quantitative 
and parameter-free agreement between simulation data 
and theory could be achieved over six time decades below 
the filament's longest relaxation time. In the "quench" 
case, short-time deviations could be attributed to the 
finite backbone extensibility of the bead-spring chains 
used in the simulations. Finally, we discussed quantita- 
tive implications for possible experimental realizations. 
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adapted a widely used scaling argument for the onset of 
the stretch-coil transition for flexible polymers to the 
semiflexible case (dominated by bending energy), and 
commented on hydrodynamic interactions. We hope that 
our thorough discussion of the non-equilibrium dynamics 
of an initially straight polymer will help to design new 
quantitative single molecule experiments and lead to a 
better understanding of more complex phenomena such 
as force transduction and recoil of disrupted stress fibers 
in cells [131 • 
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